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We studied the consequences of the Navl.4 mutation R1448H 
that is situated in the fourth voltage sensor of the channel and 
causes paramyotonia, a cold-induced myotonia followed by 
weakness. Previous work showed that the mutation uncouples 
inactivation from activation. We measured whole-cell Na^ cur- 
rents at 10, 15, 20, and 25°C using HEK293 cells stably trans- 
fected with wildtype (WT) and R1448H Na+ channels. A Markov 
model was developed the parameters of which reproduced the 
data measured on WT and R1448H channels in the whole volt- 
age and temperature range. It required an additional transient 
inactivated state and an additional closed-state inactivation tran- 
sition not previously described. The model was used to predict 
single-channel properties, free energy barriers and temperature 
dependence of rates. It allowed us to draw the following conclu- 
sions: i) open-state inactivation results from a two-step process; 
ii) the channel re-openings that cause paramyotonia originate 
from enhanced deactivation/reactivation and not from destabi- 
lized inactivation; iii) the closed-state inactivation of R1448H is 
strikingly enhanced. We assume that latter explains the episodic 
weakness following cold-induced myotonia. 

Key words: Paramyotonia, markov model, sodium channel, closed- 
state inactivation, channelopathy, skeletal muscle 



Introduction 

Paramyotonia congenita (PC) is characterized by 
muscle stiffness provoked by exposure to cold and par- 
ticularly by exercise in cold environment (1). During 
deep cooling the myotonia disappears and gives way to 
flaccid paralysis which may last several hours. Causa- 
tive mutations are in the skeletal muscle sodium channel 
Navl.4. Investigations of the biophysical alterations in 
channel gating due to PC mutations has revealed several 
gating defects consistent with membrane hyperexcitabil- 
ity. Mutant channels inactivate more slowly and with less 
voltage dependence than WT channels, deactivate more 



slowly, and exhibit a more rapid rate of recovery from fast 
inactivation (2). The very frequently occurring R1448H 
mutation which affects the outermost amino acid of the 
transmembrane segment S4 of domain DIV has been at- 
tributed to an uncoupling of fast inactivation from activa- 
tion (3). 

Voltage-gated Na+ channels are essential for the gen- 
eration of action potentials. They consist of four homolo- 
gous domains (DI to DIV) which each contain six trans- 
membrane segments (SI to S6). At depolarization, the 
S4 segments, which contain several positive amino acid 
residues and therefore function as voltage sensors, can 
move outwardly and thereby alter channel confirmation 
and function. Different charge contents of the various S4 
segments suggest that the charges have domain- specific 
functions. While S4 of DI and DII are thought to play a 
prominent role in Na^ channel activation, S4 of Dili and 
DIV regulate fast inactivation (4). Finally, the pore with 
its selectivity filter is lined by the loops between S5 and 
S6 and the S5 and S6 segments itself. 

Na+ channel activation is a multi-step process which 
is usually implemented as a series of closed states lead- 
ing to one or more open states. Generally, the distribu- 
tions of single-channel open times follow a single ex- 
ponential (5). Inactivation is coupled to activation (6). 
The voltage dependence of inactivation is supported by a 
significant gating current during fast inactivation (7) and 
can be used to describe voltage dependence of channel re- 
opening and mean open times (8). The time course of fast 
inactivation is reported to be single-exponential as well as 
double-exponential which is either implemented as two 
open states (9) or by a two-step inactivation process (10). 
After reaching the fast inactivated state, Na+ channels do 
not go immediately back to the closed states, repolari- 
zation of the membrane is necessary to initiate recovery. 
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To account for recovery from fast inactivation, which is 
not occurring by re-entering the open state (11), models 
were expanded with transitions between inactivated and 
closed states (12). As inactivation occurs from open as 
well as from closed states, and recovery from fast inac- 
tivation develops with a delay (13), multiple inactivation 
states are assumed. 

Since low temperature is the trigger for paramyo- 
tonia, temperature effects have been studied and shown 
to affect both the kinetic and steady- state parameters 
of Navl.4 WT and R1448H channels. This is not sur- 
prising, given that each of the voltage-dependent gat- 
ing steps is likely to involve different conformational 
changes in the channel and so require the breaking and/ 
or forming of chemical bonds with different energies. 
However, data obtained at room temperature cannot be 
extrapolated to physiological temperatures using a sin- 
gle temperature scaling factor. Therefore measurements 
in a wide temperature range and a suitable gating model 
which is valid in a large potential and temperature range 
are required to study R1448H. In the present study, we 
characterized the gating of Navl.4 WT and R1448H 
mutant channels with the whole-cell configuration of the 
patch-clamp technique between 5 and 30 °C. Also, we 
determined parameters of a Markov model which was 
able to fit the measurements at all potentials and tem- 
peratures. The model was then used to predict gating 
currents and single-channel properties. 

Materials and methods 

Na^ channel expression 

WT and mutant (R1448H) a-subunit constructs of 
human skeletal muscle Na^ channels were assembled in 
the mammalian expression vector pRC/C MV and trans- 
fected into human embryonic kidney cells (HEK 293) 
by the calcium phosphate precipitation method. Since 
transient expression was low (< 10%) stable cell lines 
were obtained by antibiotic selection as previously de- 
scribed (14). 

Recording techniques 

Whole-cell currents were recorded using an Axo- 
patch 200 A patch-clamp amplifier (Molecular Devices, 
USA). Signal acquisition and processing was done using 
the DigiData card (1200) and pCLAMP (V6) software 
(Molecular Devices, USA). Whole-cell currents were 
filtered at 10 kHz, and digitized at 10 or 20 jis. Patch pi- 
pettes were pulled on a Zeitz Puller (Zeitz Instruments, 
Martinsried, Germany). Pipette resistance ranged from 
0.8 to 1.2 MO. The extracellular recording solution was 
(in mM): 150 NaCl, 2 KCl, 1.5 CaCl2, 1 MgCl2 and 10 



HEPES, titrated to pH 7.4 with NaOH. The pipette solu- 
tion was (in mM): 105 CsF, 35 NaCl, 10 EGTA and 10 
HEPES, titrated to pH 7.4 with CsOH. After achieving 
the whole-cell configuration, cells were held for 10 min 
at -120 mV to ensure proper diffusion of the pipette solu- 
tion into the cell and to stabilize of the Na+ current ampli- 
tude. Recording temperature was maintained between 5 
and 30°C by a Peltier device and a HC-IOOA temperature 
controller (Dagan, USA). To avoid evaporation at high 
temperatures and dilution by condensation at low tem- 
peratures, the bath solution was continuously exchanged 
by a gravity driven perfusion system. 

Electrophysiological protocols and data analysis 

Voltage dependence of activation was obtained by a 
series of 50 ms depolarizing pulses from a holding poten- 
tial of -140 mV ranging from -85 to 55 mV and steady- 
state fast inactivation was obtained by 200 ms conditioning 
pulses from -150 to -45 mV from a holding potential of 
-140 mV followed by a test pulse to -15 mV. Activation and 
steady-state fast inactivation curves were fit with standard 
Boltzmann function as previously described (15). Time 
constants of fast inactivation were obtained by fitting dou- 
ble-exponential functions to the decaying part of the cur- 
rent traces obtained with the activation protocol. Because 
the fast component accounted for > 90% of the current 
amplitude, macroscopic inactivation of the Na^ current 
was quantified by the fast component only. Time course 
of entry into fast inactivation (closed-inactivation) was ob- 
tained by a double pulse protocol. From a holding potential 
of -140 mV a conditioning pulse V,,„d (-100, -90, -80, -70 
mV) for increasing durations (from 0. 1 to 300 ms) was ap- 
plied in order to inactivate Na+ channels without opening. 
The conditioning pulse was followed by a test pulse to -15 
mV to determine the fraction of non-inactivated channels. 
Time course of entry into fast inactivation was obtained 
by fitting a single exponential function to the normalized 
curve. Recovery from fast inactivation was determined by a 
double pulse protocol. A 150 ms pulse to -15 mV was used 
to inactivate all Na+ channels. A test pulse to -15 mV fol- 
lowed after an increasing interval (from 0.025 to 250 ms) 
at the recovery potential (-140, -120 and -100 mV). Time 
course of recovery from fast inactivation was obtained by 
fitting single/double exponential function to the normal- 
ized curve. 

Curve fits and data analysis were performed with 
pCLAMP 8.0 (Molecular Devices), Excel (Microsoft, 
Inc. Redmond, WA), and Origin (MICROCAL Software, 
Inc., Northhampton, MA). Differences from WT and mu- 
tant were considered as significant at p < 0.05 (Student's 
t-test). Grouped data are presented as mean ± SEM. SEM 
is represented in graphs as bars when in exceeds the size 
of the symbol. 



23 



Boris Holzherr et al. 



Sodium channel gating model 

Recordings from activation, steady-state fast inac- 
tivation, entry into closed- state inactivation and recov- 
ery from fast inactivation were simultaneously fit to a 
gating model using an advanced version of lonFit soft- 
ware (16). Model parameters were optimized using the 
least squares method. Ionic currents were simulated by 
solving master equations of a continuous-time Markov 
process, 



d 
di 



(1), 



r,iT) = —-expl^ J (5) 

and can be used to determine AH and AS. 

Rate constants were used to calculate single channel 
properties. If a channel opens, the number of openings 
before inactivation follows a geometric distribution (18), 
the mean of which may be calculated from the model's 
rate constants 



N = 



^ (a3 + a2+^i) ' (^6 



(6) 



whereby Pi(t) denotes the population of state / at a 
given time t, and r,y(V) denotes the rate constant for the 
transition from state / to j. The voltage dependent forward 
r^y(V) and backward ry^(V) transition rates between state / 
and j were assumed to be single-exponential functions of 
voltage (17), 



rij(V) = r[j'exp[ 



^zjrr.y • FV 



rji(y) = r'ji ' exp 



(2), 
(3) 



whereby vcr^j and zxrj^ represent the effective charge 
moving from an original state to the barrier peak, as a 
product of the total charge moved and the fraction of the 
electric field where the barrier peak was located, r/j and 
r/i represent the rate constants at 0 mV, including enthal- 
pic and entropic factors, /^represents the Faraday constant, 
R the ideal gas constant, V the membrane potential and T 
the absolute temperature. The initial state populations were 
determined as a steady-state solution of Eq. 1 at a holding 
potential V^oid with dPi(t)/dt=0. For steady-state fast inac- 
tivation curve, recovery from fast inactivation and entry 
into fast inaction, currents were simulated according to the 
pulse protocols and the respective current peak amplitudes 
were determined. Data sets used to determine model pa- 
rameters consisted of six current traces for test pulses of 
-40 to 10 mV, the steady state inactivation curve between 
-160 and -45 mV, time course of entry into fast inactivation 
at four different prepulse potentials (-100 to -70 mV) and 
time course of recovery from fast inactivation at three dif- 
ferent recovery potentials (-140 to -100 mV). 

To describe the energy profile, the rate constants in 
Eq. 2 and Eq. 3 were written with explicit entropic AS 
and enthalpic AH terms. The voltage independent parts 
are equal to the pre-f actors r/j and 



koT ( 



RT 



(4) 



The mean open time of single channels of the 
model was estimated by the reciprocal sum of the rates 
leaving the open state 



To 



1 



(7) 



To test the hypothesis of an increased probability of 
O^C4^l2 transitions, the steady- state probability was 
calculated by 



P{0 



P2 



"3 



/?2 + «6 «3+«2+/?l 



(8) 



It is very likely that there are variations in basic 
properties of channel population from cell to cell, and 
this variation may mimic the real variation seen in native 
preparations. For this reason all fits and simulations were 
done by using data of individual cells and results were 
pooled afterwards. 

Results 

Whole-cell currents 

At all temperatures activation kinetics and sodium 
currents decay were slower for R1448H than for WT 
(Fig. lA). Cooling from 30°C to 10°C slowed kinet- 
ics ~ 10-fold and reduced peak current amplitudes to a 
quarter for both WT and mutant channels (Fig. IB). In 
contrast, cooling increased the total sodium influx in- 
to the cell by different amounts: at 10°C in relation to 
30°C, the area under the curve was multiplied by a fac- 
tor of two for WT and by a factor of four for R1448H 
(Fig. IC). 

Steady-state activation curves were almost identical 
for WT and R1448H regardless of temperature (Fig. 2A, 
Table 1). Cooling decreased activation slope factor from 
~-7mV to —lOmV and potentials at half maximal ac- 
tivation were shifted by ~+8 mV to the right for WT 



24 



A gating model for wildtype and R1448H Nav1 .4 channels in paramyotonia 





r- 




Kre yrc 



IPC iFC M'C 



Figure 1- A Raw data. Representative whole-cell current 
traces recorded at different temperatures from HEK293 
cells stably expressing either WT (left) or R1448H (right) 
mutant channels: 10°C (top), 20°C (middle) and 30°C 
(bottom). Note the slowed inactivation of the mutant. B 
Temperature effect on amplitude. Temperature depend- 
ency of currents through WT or R1448H Nav1 .4 channels 
normalized to values at 30°C. C Temperature effect on 
flux. Temperature dependency of Na+ influx through WT 
or R1448H Nav1.4 channels. For B and C, values are 
mean ± SEM (n = 6). SEM is shown as bars and * indi- 
cates a significant difference between WT and R1448H 
(p < 0:05). Note that while the amplitude shows similar 
temperature dependence, the flux of the mutant is in- 
creased due to the slowed inactivation. 



and R1448H alike. Rise time of activation at 0 mV and 
higher was significantly increased in R1448H compared 
to WT (p < 0.05, Fig. 2B). Steady-state inactivation dif- 
fered significantly (p = 0.05) for the mutant as well: 
R1448H curves were significantly shifted to the left by 
~6 mV and revealed an increase of slope factor by ~4 
mV (Fig. 2 A, C, D, Table 1). Since deactivation cannot 
be measured at room temperature, we cooled to 15°C, 
10°C and 5°C to resolve sufficient data points for a fit. 
Deactivation time course was almost indistinguishable 
for mutant and WT except for the near-threshold voltage 
of -70 mV (Fig. 2E). 

For threshold-near potentials, the time constants of 
fast inactivation Tj, form the open state were smaller for 
R1448H than WT while at more depolarized potentials, 
they were larger than for WT (Fig. 3: OSI). The difference 
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Figure 2. A Activation and steady-state fast inactiva- 
tion. Activation and steady-state fast inactivation curves 
for WT and R1448H. Voltage dependence of activation 
was determined by 50 ms depolarizing pulses to the in- 
dicated potentials from a holding potential of -140 mV. 
Steady-state inactivation curve was determined with a 
200 ms prepulse to the indicated potentials prior to a 10 
ms test pulse to -15 mV (-140 mV holding potential). Solid 
lines represent fits to standard Boltzmann functions. Note 
that activation is similar but inactivation is changed for 
the mutant. B Rise time. 10-90% rise time as a measure 
for activation WT and R1448H. Note the slower kinetics of 
the mutant. C and D parameters of inactivation. Tempera- 
ture dependence of Boltzmann parameters of WT and 
R1448H for steady-state fast inactivation. Note the left- 
shift (more negative V0.5) and decreased steepness (larg- 
er k value) of the mutant. E Deactivation time constant. 
Comparison of deactivation time constants at -70 mV, the 
only potential with a significant difference. Note the ac- 
celerated kinetics of the mutant. For A to E, all values are 
mean ± SEM (n = 8-30). SEM is shown as bars when it ex- 
ceeds the size of the symbol and * indicates a significant 
difference between WT and R1448H (p < 0:05). 



in time constants was especially prominent in the voltage 
range of -60 to -30 mV and markedly increased with cool- 
ing. Cooling slowed fast inactivation of WT and R1448H 
at all voltages tested and shifted the point of intersection 
of WT and R1448H curves to more negative potentials. 
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Table 1. Boltzmann parameters of G(V) and SSFI curves. 







5 


10 


15 


20 


25 


30 


WT G{\J\ 

Tel V_/ 1 V J 




-31 .7 ± 1.6 


-35.7 ± 1.4t 




-'^8 fi + 0 7 


-40.4 ± 1 0 


-37.6 ±1.2 




k(mV) 


-9.7 ± 0.2 


-8.3 ± 0.2t 


-7.8 ± 0,1 


-7,2 ± 0.1 


-6.8 ± 0.2 






n 


S 


13 


20 


30 


19 


14 


R1448H G(V) 


Vo,{mV) 


-31.3 ±1.3 


-37.0±1.0t 


-38.3 ± 07 


-40.3 ±07 


^1.5± 1.1 


-39.4 ± 0.5 


k(nnV) 


-9.6 ±0.2 


-8.3 ± 0.2t 


-8.0 ±0,1 


-7.3 ±0.1 


-6.9 ±0.3 


-7.0 ±0.1 




n 


13 


18 


12 


19 


10 


14 


WTSSFI 


Vo5{mV) 


-88.2 ± 0 J 


-87.3 ±07t 


-84,8 ± 0.4 


-83.0 ±0.5 


-82.4 ±0.9 


-79.9 ± 1. Of 




k(mV) 


8.06 ± 0.2 


6.8 ± 0.2t 


6.4 ±0.1 


6.1 ±0.1 


6.1 ±0.1 


6.2 ±0 1 




n 


8 


13 


19 


30 


19 


14 


R1448H SSFI 


Vo5{mV) 


-93.7 ± i.r 


-91.0±0.8*t 


-887 ± 1.2* 


-87.0 ± I.r 


-867± 1.1* 


-85.0 ± 1.3* 




k(mV) 


12.0 ±0.3* 


10.7±0.2*t 


108±03* 


10.1 ±0.2' 


9.9 ±0.4" 


9.5 ± 0.3* 




n 


13 


18 


12 


18 


10 


14 



* indicates a significant difference between WT and R1448H at tlie same temperature, f indicates a significant difference between 10, 
20 and 30°C, (p < 0:05). 
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Figure 3. Tinne constants. Time constants from and into 
the fast inactivated-state were plotted against the corre- 
sponding membrane potentials. Recovery, entry (Closed- 
state inactivation, CSI) and inactivation from the open- 
state (OSI) were determined forWT and R1448H between 
5 and 30°C. Values are mean ± SEM (n = 5-30). SEM is 
shown as bars when it exceeds the size of the symbol 
and * indicates a significant difference between WT and 
R1448H (p < 0:05). Note the accelerated CSI and slowed 
OSI of the mutant at all temperatures, while the recovery 
is indistinguishable. 



Additionally, R1448H reduced voltage dependence of 
for all temperatures tested. 

R1448H accelerated entry into closed- state inactiva- 
tion (CSI) by about two-fold on average (Fig. 3: CSI, Ta- 
ble 2). The left-shift of the steady-state inactivation curve 
may explain this enhanced closed-state inactivation. The 
mutation reduced its voltage dependence, possibly by the 
removed S4 charge, and slowed the open-state inactiva- 



tion. Finally, R1448H showed a tendency to recover more 
rapidly without reaching significant levels (Fig. 3: recov- 
ery. Table 3). However, neither the delay in onset to re- 
covery nor its voltage and temperature dependence were 
altered by R1448H. 

Gating model 

The gating model used in the present study consisted 
of a series of four closed states C1-C4, one open state O 
and four inactivated states I1-I3 and Ij (Fig. 4). By conven- 
tion, all transitions towards O have positive valences be- 
cause they are favored by depolarization, while those away 
from O have negative valences because they are favored 
by repolarization. Rate constants between the closed states 
and rate constants between C3-C4-O and Irl2-l3 were as- 
sumed to be equal to reduce the number of free parameters. 
The model is based on previous models of Vandenberg and 
Bezanilla (19). However, to properly describe gating of 
WT and R1448H and their temperature dependence, two 
modifications were made. First, we introduced a transient 
inactivated state Ij as previously suggested by Kahlig et 
al. (10) to account for the biphasic inactivation especially 
at low temperatures. Second, we introduced a transition be- 
tween C4 and I2. This transition was essential to reproduce 
inactivation from closed states especially in the voltage 
range around threshold of Na+ channels. 

The model was able to reproduce all measurements 
including the strong voltage dependence of channel acti- 
vation, open and closed state inactivation, recovery and 
temperature dependence (Fig. 5). The model resulted in 
rate constants of WT and R1448H which were similar 
for most transitions (Tables 4 and 5). However the rate 
constants alpha3, beta3 and alpha6 markedly differed 
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Table 2. Time constants of entry into fast inactivation. 







5 


10 


15 


20 


25 


30 


WT 


^en(ry-10o('^^) 


47.0 ±7.4 


32.5 + 3.8t 


16.1 +1,6 


10.1 +1,3 


5,9+0.8 


3.4+0.7t 




n 


8 


15 


18 


22 


15 


8 




'^eillry-9o('^S) 


70 8 ±4.8 


51 3±2.7t 


31 9 ± 1.9 


19 7 ±0 .9 


10 3 ±0.5 


5.2 ±0.7t 




n 


7 


15 


20 


24 


15 


9 




'^eri<fy-eo('^^) 


58 7 + 4 0 


40.4 + 2.4t 


32 9 + 2 1 


23 5±09 


12 9 ±0.7 


S O ± 0.5t 




n 


7 


15 


20 


24 


15 


9 




'^an(iY-7o(ms) 


38.1 ±2.6 


26.2±1.8t 


22.1 ±1.4 


17.1 ±1.0 


10.0 ±0.8 


7.6±0.8t 




n 


7 


13 


20 


22 


14 


9 


R1448H 




264+1.2* 


14 1 +0.7*t 


9.2+0.5* 


5.5+0.3* 


3.2 ±0.2* 


2.0 ± 0.2t 




n 


10 


18 


12 


17 


10 


9 




^9nliY-9o('^S) 


27.8 ± 0.9* 


17.0±0.7*t 


11.2 ±0.5* 


6.9 ±0.3* 


4.3 ±0.2* 


2.6 ±o.rt 




n 


10 


18 


12 


16 


9 


9 






25.3 ±0.7* 


15.8±0.6*t 


10.9 ±0,5* 


7,2 ±0,2* 


4.4 +o.r 


2.9±0.2*t 




n 


10 


18 


12 


17 


9 


9 






20 1 ±0.7* 


12 3±0.5*t 


8.7 ±0.4* 


6.0 ±0.2* 


3.5 ±0.1* 


2.6 ±0.2*t 




n 


10 


18 


12 


17 


8 


9 



* indicates a significant difference between WT and R1448H at the same temperature, f indicates a significant difference between 10, 
20 and 30°C, (p < 0:05). 



Table 3. Time constants of recovery from fast inactivation. 



X 




5 


10 


15 


20 


25 


30 


WTf140 mV) 


delay (ms) 
n 


16,61 ±0.49 

1.54 ±0.03 
6 


7,94±0.42t 

0.81 ±0.04t 
14 


3.28 ±0.11 
34.6 ± 3.8 
0.57 ±0.01 

17 


1.37 ±0,04 
24,5 ± 2,2 
0.31 ± 0.01 
19 


0,66 ± 0.03 
19,1 ± 1.4 
0,18 ±0.01 
10 


0,35 ± 0,03t 
16,6 ±0,8 
0.1 ±0.003t 
9 


R1448H('140 mV) 


j(ms) 
Tret ^(ms) 
delay (ms) 
n 


15.50 ± 1.79 

1.79 ±0.25 
9 


6.45 ± 0.28t 

0.84 ± 0 03t 
17 


3.23 ±0.14 
39.3 ± 4.1 
0.58 ± 0.02 
12 


1.55 ±0.08 
32.6 ± 3.6 
03±0.01 
18 


0.76 ± 0.04 
21,1 ± 1.3 
0.19 + 0.02 
7 


0.45 ± 0.04t 
21 5 ±1.9 
0.099 ± 0.003t 
7 


WT(-120mV) 


Tret j(ms) 

Trec_z(ms) 

delay (ms) 
n 


40 92 ± 1.09 

3.26 ± 0.07 
6 


19.14 ± 1.16t 

1.51 ±0 lot 
14 


7.81 ± 0.36 
63.1 ±6.7 
1.03 ± 0 02 
18 


3 21 ± 0.11 
50,2 ± 5.6 
0 55 ± 0.02 
19 


1.45 + 0.09 
38,2 ± 3.2 
0 29 ± 0.01 
10 


0.75 ± 0.06t 
36,0 ±3.0 
0.15 ±0.01t 
9 


R1448H(-120 mV) 


Tret_i(ms) 

Tr^_2(ms) 

delay (ms) 
n 


31.44 ±3.26* 

3.34 ±0,48 

9 


13.86 + 0 66-t 6.88+0.28 
79.5 ±15 9 
1.45±0.05t 0.97 + 0 04 
17 10 


3 18 ± 0.18 
51.&±6.S 
0,57 ±0.02 
18 


1 60 + 0.16 
39 4 ± 3.4 
0,30 ± 0.01 
6 


0.95 ± 0.07't 
47 8 + 3.5 

0,14±0,01t 
7 


WT(-100 mV) 


Tret (nns) 

Tr«_2(mS) 

delay (ms) 
n 


59.97 ± 2.89 

4.20 ±0.15 
5 


52.82 ± 2.90t 

3.31 ±0 23t 
13 


26.21 ±1.40 

1.93 ±0 07 
19 


9 08 ± 0.35 
58,2 ± 5,8 
1 10 ± 0.04 
16 


4.72 ± 0.34 
98,6 ± 26,2 
0.56 ± 0.03 
8 


2.27 ± 0.24t 
96,0 ±18,5 
0.33 ± 0.04t 
10 


R1448H(-100 mV) 


Tf« i(nns) 
tr« 2(ms) 
delay (ms) 

n 


59.30 ± 3.32 

6.01 ±0 78 
8 


29.17 ± 1.43't 16.85 ±0 86* 

3.14 ±0 19T 1 79 ± 0 14 
17 12 


8 36 ± 0.43 

0 96 + 0.06 
17 


3.75 ± 0.30 
53 2 ± 4.4 
0 58 ± 0.05 
6 


2.10 ±0.1 1t 
110.8 ± 19.5 
0.27 ± 0.03t 
8 



* indicates a significant difference between WT and R1 448H at the same temperature, f indicates a significant difference between 1 0, 
20 and 30°C, (p < 0:05). 



between WT and R1448H. The smaller alpha6 explains alphaS and betaS, the transition between C4 and 12. In 
the impaired fast inactivation of R1448H. The most contrast to the WT for which the calculated rate con- 
striking difference between WT and mutant related to stants suggest that this transition does not occur, the mu- 
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tant performed this closed-state inactivation transition 
with a high likelihood. 

Effective charge movement increased with tempera- 
ture for all transitions in both WT and R1448H by about 




Figure 4. Gating model. Gating model used for kinet- 
ic simulations, consists of 4 closed-states 01-04, one 
open- state O, four inactivated-states 11-13 and IT. Ar- 
rows between states indicate possible transitions and J 
and J (i=1 ..6) represent the voltage-dependent rate con- 
stants. Rate constants _3 and _4 were calculated based 
on microscopic reversibility. Note the additional transient 
inactivated state IT and the transition 04-12 both of which 
have not been previously used. 



30%. For example the equivalent gating charge for open- 
ing the channels from rest was 5.31 to 6.78 eO for WT 
and 4.39 to 5.56 eO for mutant, obtained by summing the 
charges for each of the four transitions from CI to O (val- 
ues dependent on voltage, Table S6). Most of the charge 
movement in the activation pathway was concentrated 
in the last transition (C4-0) 2.58 ± 0.06 to 3.06 ± 0.04 
eO for WT and 2.53 ± 0.05 to 2.98 ± 0.09 eO for mu- 
tant. We interpret this finding so that this transition may 
represent several steps in one the final of which may re- 
ally be voltage-independent. In general, effective charge 
movement for transitions from inactivated to closed states 
during recovery were notably larger compared to their re- 
spective forward rates during closed- state inactivation ac- 
counting for the strong voltage dependence of recovery 
from inactivation. For the mutant the equivalent gating 
charge movement during recovery was smaller than for 
WT leading to reduced voltage dependence. About 50% 
of total gating charge of WT and 40% of charge for the 
mutant was immobilized by fast inactivation. 

Free energy barriers 

The energy changes involved in the transitions be- 
tween the closed-states (C1-C2-C3-C4) and the paral- 



Table 4. Model parameters for WT 



Transition 


Parameter 


10X 


IS'C 


20"C 


25X 


C - C and li - b 


al 


0.1957 ±0.0020 
7.75e03±3 10e02 
-0 7134 + 0.0144 
1.61e03±6 90e01 


0.183210.0132 
1.25e04±1.96e02 
-0.8312 + 0.0205 
1.See03±1.09e02 


0.171610.0153 
2.27eO4±3.30e02 
-0.9831+0.0166 
1.87e03±1.09e02 


0.167810.0125 
3.36e04±1 iae03 
-1.0712 + 0.0167 
2.19e03±6 99e0l 


C4- 0 and i^- I3 


zx a2 
ot2 

^2 


1,907910.0422 
2.eOe04+1 46e03 
-0.6711+0.0374 
5.93e01±5.17e00 


1.964910,0323 
6 37e04±3.65e03 
-0.8309 + 0 0596 
1.See02±4.02e01 


1.977510.0164 
1.43e05 + fl4ie03 
-0.9167 + 0.0307 
570e02±2.44e01 


1,897810.0341 
3.50e05+1.30e04 
-1.1666 + 0.02B9 
2.03e03l3.29e01 




zxaa 
03 

im 

n3 


0.0001 ±0.0001 
5.40e-07 + 1.61e'07 
*1 .098210.0034 
3.47e-10±1.09e-10 


0.015310,0122 
5 56e-06±1.53e-06 
•1.219010.0254 
3.816-091 1.1 6e-09 


0,0552 + 0.0475 
3.l4e-05±1 57e-05 
■1.343710.0406 
2.22e-0e±1.35e-O8 


0.031010,0201 
5.14e-04±1.39e'04 
•1,5462 ±0.0251 
2.74e-07±6.99e-03 




zxa4 
04 

P4 


0.000410.0002 
5.36e01±5.63e00 
-0.846410,0113 
e.56e-01 + 1.76e'01 


0014010.0123 
1 19e02±8.12e00 
0.632610.0349 
3 30e00±4.11e-01 


0.020710.0160 
2.4ae02±1.06e01 
-0,953410.0291 
1.45e01±2 02e00 


0.000410.0002 
5.46e02±1 71e01 
-1,073110,0304 
4.79e01±4 51eO0 




zxaS 
05 


0.0016 10.0006 
3.78e02i1.76e01 
-1.096710.0032 
2.40e-01 ± 1 70e-02 


0.004010.0026 
4.74e02±2.30e01 

-1 2302 + 0.0258 
3 06e-01±3.96e-02 


0.011310.0076 
4.92^02 ±2,26e01 
-1.387710.0310 
3.19e-01±4 81e-O2 


0.012310.0107 
6.89*02 t2,06e01 
-1,565010.0088 
3.73e-01 ± 1 60e'02 


O-It 


2xa6 
06 

P6 


0.251410.0098 
5.29e02±2 24e01 
-0.0002 + 0.0001 
2.40e01±2 24e00 


0.334910.0148 
1.42e03 + 5.73e01 
-0.002710.0016 
3 27e01±1.47e00 


0.406310.0258 
3.0SeO3+1.06e02 
-0.0185 + 0.0164 
3.37eOl±1 40e00 


0.469210,0131 
6.05e03 + 3 12e01 
-0 0345 + 0.0164 
3.90e01±2 65e00 



Model parameters ai and (3 i (i = 1 ..6) are the corresponding forward and backward transition rate constants at 0 mV. Parameters zxai 
and zx(3i are the valences of the corresponding transitions and reflect the voltage dependence of transition rates.* indicates a signifi- 
cant difference between WT and mutant (p < 0:05). 
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Figure 5. Representative current traces and fits. Original data at 10 to 25°C (black dots) are superinnposed with fits of 
the model (red lines) obtained by simultaneous fitting of A) a series of 6 current traces elicited by pulses from -40 to 10 
mV, (B) the steady-state fast inactivation curve determined in the range between -150 and -40 mV, (C) the time course 
of entry into fast inactivation at prepulse potentials of -1 00, -90, -80 -70 mV and (D) the time course of recovery from fast 
inactivation at recovery potentials of -140, -120, -100 mVfor WT (top) and R1448H (bottom). Note the quality of the fits. 
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Table 5. Model parameters for R1448H. 

Transition Parameter 10'C 



20X 



25X 



C - C and li - Ig 



zxal 
al 



0.133910.0145* 
7.18e03±1.61e02 
-0.4894 ± 0.0168- 
3.06e03±1.14e02* 



0.1122±0.0083* 
1.27eO4±4.12e02 
-0.6171 ±0.0192* 
4.06e03±1,77e02* 



0.126410.0129* 
2.26e04±7.90e02 
-0.6763 ± 0.0177* 
4.97e03±9.29e01* 



0.116610.0165* 
3.62e04±1.60e03 
-0.7396 10.0245* 
6.098031 5.45e02* 



C4-Oandl2-l3 zx ci2 1.858410.0441 1.897010.0513 1.900110.0469 1.804110.0828 

a2 2.93e04l1.60e03 7.08e04l6.47e03 1.52e05l7.35e03 4.27e05l 3.24e04* 
zx32 -0.6685 10.0310 -0.7235 10.0472 -0.8106 10.0200* -1.1 766 10.0449 
P2 2.46e01l2.65e00* 1.09e02l5.04e00* 3.73e02t 2.25e0r 1.06e03l 1.70e02* 

C4-I2 aa3 0.002710.0010 0.0291 10.0278 0.000910.0003 0.003610.0020 

a3 7.50e02±3.55e0r 8.37e02l 5.29eOr 8.07e02±4.37e01* 1.07e03l4.22e01* 

zx^S -0.9188 ±0.0187* -0.9477 ± 0.0324* -1 .0421 ± 0.0270* -1 .2099 ± 0.0254* 

1.19e00±1,17e-0r 1 .28eO0± 1 .65e-0r 1.39e00± 1.69e-01* 1 .92e00l2.54e-0r 



It- I3 


zxa4 


0.0051 ±0.0024 


0.0022 ±0.0012 


0.0135±0.0124 


0.0009 ±0.0003 




a4 


3.53e01 ±3.88e00* 


5.53e01±2.73e00* 


1.39e02±1.91e0r 


1.87e02±8.12e00* 




zxp4 


-0,6289 ±0.01 76* 


-0.6326 ±0.0283* 


-0.6043 ±0.0322* 


-0.7623 ±0.0335* 




,■4 


4.286-01 ±5.99e-02- 


1.12eO0±1.35e-0r 


5.9l600± 1.21600* 


1.38601 ±2.93600* 


C3-I1 


zxu5 


0.0036 ±0.0024 


0.0020 ±0.0011 


0.0001 ±0.0001 


0.0013 ±0.0008 




a5 


5.34e02±7.97e01 


7.25e02±6.43e0r 


8.75602 ±3.03e0r 


1. 08603 ±5.57e0r 




zx^S 


-0.917910.0183* 


-0.9749 ±0.0147* 


-1.04291 0.0270* 


-1.2121 ±0.0247* 






8.03e-01±8.72e-02* 


1.06e00±8.29e-02* 


1.52e00±1.94e-0r 


1.89e00±1.89e-0r 


0 - f. 


Z>:i(6 


0 2861 ±0.3071- 


0.332810.0066 


3 425010.0110 


0 4500±0.0129 




u6 


1.16e02±7.78e00- 


2.61e02± 1.04e0r 


5.56e02± 1.47e0r 


1.03e03± 3.54601* 




zx^6 


-0.001 5 ±0.0012 


0.0093 ±0.0092 


-0.0002 ±0.0001 


-0.0001 ±0.0001 




P6 


1.61e01±2.11e00* 


2.0 1e01 ±1.54600* 


2.39601 ±2.96600* 


2.75601 ±3.06600* 



r 8 8 3 6 



Model parameters al and (3 i (i = 1 ..6) are the corresponding forward and backward transition rate constants at 0 mV. Parameters zxal 
and zx(3i are the valences of the corresponding transitions and reflect the voltage dependence of transition rates.* indicates a signifi- 
cant difference between WT and mutant (p < 0:05). 



Table 6. Equivalent gating charges. 

WT Transitions No 10 15 *C 20 °C 25X 



C - C and Ih| — I2 


1 


0.91 ±0,01 


1.02 ±0,02 


1,15±0.02 


1.24 ±0.02 


C4 - 0 and I2 - h 


2 


2.58 ±0.06 


2.80 ± 0,07 


2,69 ±0,03 


3.06 ±0.04 




3 


1.10±0.01 


1.23± 0 03 


1.40 ±0 06 


1.58 ±0.03 




4 


0,85 ±0.01 


0.90 ± 0,04 


0 97 ±0 03 


1,07 ±0.03 


C3-I1 


5 


1.10 + 0.01 


1.23± 0 03 


1.40 ±0.03 


1.58 ±0.01 


0-It 


6 


0.25 ±0.01 


0.34 ±0.01 


0,42 ±0.03 


0.50 ±0.02 


R1448H Transitions 


No 


10 X 


15 X 


20 X 


25X 



C - Cand I, - 
C4 - O and \2 
C4-I2 

C3-I1 

o-h 



1 
2 
3 
4 
5 
6 



0,52 ±0.02 073 ±0.02 0.80 ±0.02 0.86 ±0.03 



2.53 ±0.05 
0.92 ± 0.02 
0.63 ± 0.02 
0.92 ± 0.02 
0.29 ± 0.01 



2.62 ±0.07 
0.98 ±0 04 
0.63 ±0 03 
0.98 ±0 01 
0.34 ±0 01 



2.71 ± 0.05 
1.04 ±0.03 

0 .62 ± 0 .03 
1.04 ± 0.03 
0.43 ± 0.01 



2.98 ±0.09 
1.21 ± 0.03 

0.76 ± 0.03 
1.21 ± 0.02 
0.45 ± 0.01 



Equivalent gating charges calculated for model C14. Values are mean ± SEM (n = 6-9). 



lei inactivated- states (11-12) consist of both entropic and 
enthalpic changes, suggesting that chemical bonds are 
reforming and conformational changes of the channel are 
taking place. For the C4-0 transition there is a net de- 
crease in enthalpy along with a net decrease in entropy 
when the channel goes from the last closed state C4 to the 



open state O (Table 7). This result suggests that the open- 
ing step corresponds to a reorganization of the channel 
with a decrease in the degrees of freedom of the molecule 
giving a more ordered system in the open state. While 
the energy barrier for O-I^ was increased by 5% in the 
mutant (Fig. 6, left), the one for C4-I2 was reduced down 
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A gating model for wildtype and R1448H Nav1 .4 channels in paramyotonia 
AH AS TAS20 AGIO AG25 



c - C/l^ - 12 


WT 


56 ±6 


211 + 18 


62 


-4 + 8 


-7 + 8 




R1448H 


43 ± 3 


160 ± 11 


47 


-2 ± 4 


-4 ± 5 


C4 - 0/I2 - 13 


WT 


-52± 10 


-131 +27 


-38 


-15± 13 


-13±13 




R1448H 


-34 ±15 


-66 ± 39 


-19 


-16± 19 


-15± 19 


C3-(i 


WT 


8 ± 7 


87 + 23 


26 


-17+10 


-19 ± 10 




R1448H 


-8 ±4 


27 ± 12 


8 


-15 ± 5 


-16 ± 5 


C4 - 12 


WT 


42 ± 75 


203 ± 131 


60 


-16 + 84 


-19 ±85 




R1448H 


-6± 10 


31 ±31 


9 


-15± 13 


-16±13 




WT 


88 ±8 


336 ± 25 


99 


-7 ± 11 


-13 ± 11 




R1448H 


74 + 5 


277 + 9 


81 


-5 + 6 


-9 + 6 


h-h 


WT 


-76 + 9 


-235 + 24 


-69 


-10± 12 


-6 + 12 




R1448H 


-81 ±29 


-251 ±81 


-73 


-10 ±37 


-7 ±38 



AH (kJ /mol) and AS (J /Kmol) values were obtained by fitting Eq. 4 and Eq. 5 to a'i and (3'i (i=1 ..6) values. TAS20 (kJ /mol) was calculated 
for 20°C. AGIO and AG25 (kJ /mol) were calculated using AG = AH - TAS for 10°C and 25°C respectively Values are fit values ± SEM. 



to 50%, 50 vs. 95 kJ/mol, confirming the facilitated tran- 
sitions between C4 and 12 due to strikingly increased al- 
phas, i.e. meaning enhanced closed-state inactivation for 
R1448H (Fig. 6, right). 

Single-channel behavior 

Our finding that entry into rapid inactivation of 
R1448H was faster than for WT at threshold-near po- 
tentials (Fig. 3) was interpreted as tendency of R1448H 
channels to deactivate and inactivate through closed 
states. To further prove this hypothesis, the probability 
of transitions from O to 12 was modeled and it is obvi- 
ous that this transition occurs in R1448H and not in WT 
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Figure 6. Free energy barriers between states. Total free energy barriers be- 
tween states were calculated for -160 mV (solid line) and + 50 mV (scattered line) 
for WT (black) and R1448H (red). The value to the left of the energy barrier was 
set to 0 to allow direct comparison between curves. Note the very low barrier for 
the C4-I2 closed-sate inactivation transition (right) for the mutant. In contrast the 
0-IT open-state inactivation transition barrier is only very mildly elevated for the 
mutant (left). 



(Fig. 7). Cooling shows a clear increase in the probability 
for this transition as expected from the whole-cell current 
data at lower temperatures. 

The model's rate constants were used to calculate 
single-channel properties to determine whether the 
slowing of the current decay observed for R1448H can 
arise from longer open times or an increased number 
of openings. The estimated mean open times were up 
to 4-fold longer for R1448H than for WT. Cooling in- 
creased the mean open time of both R1448H and WT 
channels (Fig. 8 top). The bell-shaped curves showed 
open-time maxima between -50 and 0 mV. To the left 
of the maximum, the mean open time was dominated 
by the rate beta2 and to the right 
of the maximum by alpha6. This 
means that Na+ channels open 
and close several times before 
they finally enter the inactivated 
state. Importantly the calculated 
number of openings was -20% 
greater for R1448H than for 
WT (Fig. 8 bottom). Cooling 
reduced the number of re-open- 
ings for both WT and R1448H. 
In summary the slowed decay 
of whole-cell currents (Fig. 1) is 
due to an increase in open times 
which are further increased 
by cooling. The rate constants 
and the transition probabilities 
showed that the increased num- 
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Figure 7. Voltage dependence of closed-sate inactiva- 
tion probability. The probability for a transition from O 
C4 12 was calculated according Eq. 8 for WT (open 
symbols) and R1448H (filled symbols) for 10°C-25°C. 
Values are mean ± SEM (n = 6-9). SEM is shown as bars 
when it exceeds the size of the symbol. Note that this 
transition practically does not occur in WT (P~0), but only 
in the mutant. 
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Figure 8. Temperature and voltage dependence of open times and number of 
openings. Temperature and voltage dependence of the mean open time (top) 
and the number of openings before inactivation (bottom) was calculated for in- 
dicated voltages for WT (left) and R1448H (right) between 10 and 25°C. Values 
are mean ± SEM (n = 6-9). SEM is shown as bars when it exceeds the size of the 
symbol. Note the elongated open time and the elevated number of openings in the 
mutant compared to WT. 



ber of R1448H openings is due to re-openings from the 
closed state C4 and not from the inactivated states. As 
the mutant channel shows the minimum of the energy 
landscape for 13, the channels reach this state by the C4- 
12 pathway instead of by IT. Mutant channels go along 
the O^C4^I2^I3 pathway. 



Discussion 

Our whole-cell data confirms previous studies in so far 
as R1448H slows open-state inactivation and shifts steady- 
state rapid inactivation to more negative potentials (3, 20- 
22) and that the seemingly temperature sensitivity in par- 
amyotonia is a result of channel kinetics which are already 
slowed in the warmth and undergo a normal slowing with 
cooling (23, 24). Therefore, we assume that the required 
changes made to our model to best fit the data are not the 
result of our specific measurement or our set-up but rather 
reveal generally valid states and transitions. 

The required introduction of the transient inactivat- 
ed state IT into our model suggests that open-state in- 
activation may result from a two-step process. The two 
inactivation phases become more obvious at low tem- 
peratures whereas they cannot be temporally resolved at 
higher temperatures. A biphasic 
inactivation process is actually in 
agreement with the classical HH 
model and with previous single 
channel measurements (3). We 
interpret the two phases to be 
linked to deactivation and inac- 
tivation. 

The required enabling of the 
C4 - 12 transition (which really 
occurs in the mutant only and 
not in WT) is mainly responsi- 
ble for the channel re-openings 
of R1448H which are known to 
cause repetitive action potentials 
and paramyotonia. The transition 
rates and energy barriers of our 
model suggest that the re-open- 
ings originate from C4-0 transi- 
tions and not from O-I transitions 
since the inactivation on-rate is 
reduced. This view is further evi- 
denced by the enhanced inactiva- 
tion from closed states, a slightly 
accelerated recovery from rapid 
inactivation, and the absence of 
a persistent current due to the 
limited number of re-openings by 
the increased rate of C4-I2 transi- 
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tions. The enhanced deactivation has been previously also 
deduced (3). 

As found previously for R1448H but not R1448C (3), 
closed-state inactivation (CSI) is strikingly enhanced for 
the R1448H mutation. We assume this is due to a more 
outward positioned resting-state S4 because of the elimi- 
nated positive charge at residue 1448 similar to calcium 
channel mutations (16). The enhanced CSI can explain 
the transition from myotonia to flaccid muscle weakness. 
Since R1448H impairs the movement of the voltage sen- 
sor, the receptor for the inactivation gate is more read- 
ily available for voltages around the activation threshold 
and less available for further depolarized voltages. The 
slowing of the rapid inactivation prolongs the duration of 
muscle action potentials as measured in vivo (25), where- 
by the combination of repetitive activity and prolonged 
duration of each action potential leads to a cold-induced 
depolarization and thus intracellular Na+ accumulation 
that can even be detected by ^^Na magnetic resonance to- 
mography in vivo (26, 27). 
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